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1 Introduction 

Complex numbers create some of the most beautiful pictures in mathematics. 
Mandelbrot sets, Julia sets, and many other computer- generated images have 
become everyday images; hundreds of books have been published and songs 
have been written (my favorite is [9]). However, the algorithms which drive 
these computations make a major assumption about the complex plane which 
can create some visual havoc. Robert Corless ([3j) asked if there is a way to 
compensate for this error. Unfortunately, we show this is not possible. 

2 Complex Multiplication and Exponentiation 

The multiplication of two complex numbers is rather straight-forward: for 
zi,Z2 G C, where zi = a + bi and Z2 = c + di with a,b,c,d GM. we have 

Z1Z2 = (a + bi){c + di) = {ac — bd) + {ad + bc)i. (2.1) 

However, this is not computationally efficient for exponentiation; the com- 
putation of (2 -|- 3i)^ is not as quick using this definition: 

(2 + 3if = {2 + 3i){2 + 3i){2 + 3i) = ((4 - 9) + {6 + 6)i)(2 + 3i) 

= {-5 + 12i){2 + 3i) 
= {-W - 36) + {-15 + 24)1 
= -46 + 9i. 

If we convert 2-\-3i to polar form z = re*^ , we can compute this exponen- 
tial much faster. The conversion equations between rectangular coordinates 
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z = a + bi and polar form z = re*^ are: 

a = r cos{e) h = rsm{e) (2.2) 

r = ^/a? + 62 tan(e) = - (2.3) 

a 

So we convert 2 + 3i to polar form: r = ^/2'^ + = -v/l3 ~ 3.60555127. . . 
and tan(0) = 1.5 =^ ^ ?a 0.9828 or 56.31°. Now we can use some simple 
properties of exponents to get: 

(^re'^)" = rV"^ fornGN (2.4) 

So we have 

(2 + 2,if = (^e°-9828i)' = (V13)' e3{0.9828)^ ^ 46.87216658e2-94838ii7^ 

Now convert this back into rectangular coordinates: 

a = 46.87216658 cos (2.94838117) = -46 
b = 46.87216658 sin(2.94838117) = 9 

so we have (2 + 3i)^ = —46 + 9i. Perhaps this polar form conversion seems 
more tedious, but consider how much easier it makes computations with 
higher powers, like (2 + Si)^^ (which, in case you're wondering, equals 
-13, 833, 225, 534, 613, 558 - 3, 190, 610, 873, 034, 597«) . 

A closer examination of the complex exponential function reveals an in- 
teresting property: is not a one-to-one function! Using Equations 12.21 and 
12.31 we see that the periodic nature of the trigonometric functions cause there 
to be multiple (co-terminal) angles which satisfy the formulae. For example, 
1 + i can be written as v2 e** and as v2 e~*. The angles (referred to as 
the argument of the complex number; arg(2;)) differ by a multiple of 27r, and 
any number of the form v^e^^'^'^^^)* with /c G Z is also a polar form of 1 + z. 
This means that the inverse of the function is not a function, but we can 
work around this. 

The standard definition of the logarithm of z requires that we specify an 
interval of length 27r for the angles. This is called a branch of the logarithm. 
The generalized logarithmic function is defined by 

log(z) = log(re*^) = ln(r) + iB (2.5) 

Now we have a point where this function is undefined: when z = we would 
need to compute ln(0) and to find 0, which is not well defined for 0. Any 
subset of the complex plane on which we can define a (one-to-one) function 
which acts as an inverse of is called a branch of the logarithm (see any 
introductory text of complex analysis for a more detailed exposition; like 
[2], pages 38-40). This requires that there is a part of the plane "removed" 
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from the domain of log(2:); this is cahed the branch cut. It is a curve from 
z = extending out to infinity; most often it is a straight hne. When we us 
the principal branch of the logarithm, the cut is the non-positive real axis. 
This means that the angles used for 9 are in the interval (— vr, vr]. This is the 
default branch in almost every computer algebra system in the world; and is 
called the principal branch of the logarithm. This means that log(2;) is not 
analytic on the whole plane. 

When we define = e"'°s(^) for n G N we have a problem: how can 
the composition of nonanalytic functions become analytic? It has to do with 
the behavior of the branch cut. We will follow the computation through the 
following chain of compositions: 

z ^ log(z) ^ nlog(z) ^ e"^°§(^) = (2.6) 

When we take a small disk A = {z : \z — r\ < e} about a point r < in the 
negative half of the real axis, it is divided into two halves by the axis (see 
Figure [1] on page [3]) : 

the upper half-disk Ajj = Af] {z\ lm{z) > 0} 
the lower half-disk = An{z\ lm{z) < 0} 

Im 




Figure 1: A disk on the negative real axis 

These half-disks are transformed by the map log(z) into regions about 
the lines Im(z) = vr and Im(z) = — tt. The images are almost semi-circular 
regions, with center at log(r) and radius of e (see Figure [2] on page|3|). 

Now we have illustrated the first step in the chain in 12.61 Our next step 
is to multiply these points by n S N. This expands the entire picture; the 
horizontal edges of the half disks are moved to Im(2;) = invr, and the radius 
will increase by a factor of n as well. So we have the image as seen in Figure 
[3] on pageHl Our final step is to exponentiate; the result will depend on the 
parity of n. If n is odd, the result is a disk of radius centered around the 
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lm(z) = ix 
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Figure 2: The image of the disk in Figure [T] under the map log(2;) 



point on the negative real axis with magnitude r" (i.e., the center is — r"); if 
n is even, the result is a disk of radius e" centered around the point r" (the 
point on the positive real axis with magnitude r"). 



Im(z) = nix 



lm(i) = -nix 
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Figure 3: The image of the regions in Figure [2] multiplied by 



n 



What is critical in this sequence of maps is the behavior of the half 
disks in Figure [3] when acted on by the exponential. When n is an integer, 
these half disks are aligned after the exponential map: the arguments of the 
horizontal segments of the boundary are nvr and — nvr which are mapped to 
points on the real axis (see Figure U] on page [5]). When any other branch 
cut is used, the same process can be applied. The result is that the function 
defined by 

z = 



(2.7) 



is analytic and well defined for n G N. 

A technical aside: we can consider the infinite family of disks, with centers 
on the lines Im(z) = {2k + 1)tt, k £ Z, as the logarithmic image of the single 
disk in Figure [TJ If we follow this line of reasoning, the expansion by a 
factor of n is the critical step; all of the disks now have centers on lines with 
Im{z) = nkn, k G Z, and the exponential will map this family of disks (in 
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If n is odd 



If n is even 



Figure 4: The image of the regions in Figure [3] under the map 



an infinite-to-one fashion) onto a single disk as seen in Figure H] on page O 

3 Non-integer Exponents 

What happens when we try to define z"' if a is not an integer? We have the 
same sequence of composition: 



z I— )• log(z) I—)- alog(2;) I—)- e 



alog{z) 



(3.1) 



but what breaks down? 

Let us take a simple example: z^ for n G N. The process we used for the 
integer exponents can be paralleled, but we must consider the infinite family 
of disks mentioned at the end of Section [2l When we multiply by a = ^ in 
the second step, we will compress n of these disks to lie between the lines 
Im(z) = —IT and Im(z) — vr. An example (with a = 2) is shown in Figure [5] 
on page[6l This results in two distinct square root, and n distinct n^^ roots 
in general, for any non-zero complex number. Note that if n is even, then 
none of the family of disks will map onto the lines Im(z) = ibvr, and there 
will be no problem when we exponentiate. However if n is odd, then there 
are images of the branch cut (the lines Im(z) = ibnvr) which are compressed 
onto the lines lm(z) = zinr. This causes members of the family of disks to lie 
on the branch cut, exactly as in Figure [2j Thus when we exponentiate, these 
half-disks realign, and the resulting function is analytic (away from z = 0). 

However, the value a = ^ is a very special case. What about the more 
general case? The first step here works exactly as in Figure [2] on page [H a 
disk about the negative real axis is split into two half-disks and mapped to 
two half-disks about the lines Im(z) = zbvr. These regions are expanded by 
a factor of a, so the lines in figure [3] on page 3] would be labelled airi and 
—ani. When we restrict ourselves to the images contained between the lines 
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Figure 5: The images of the disk in Figure [T] under the map i log( 



Im(z) = zbvr, we have something similar to that seen in figure El on page[6l 
Notice that the disks that intersect the branch cut hues are not bisected, but 
are divided into nonequal portions; this is a critical difference. 



Im(z) = ix 
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Figure 6: The images of the disk in Figure [T] under the map alog(z) 

Now we exponentiate, and here is the rub: The images of the rays 
lm{z) = ibavr do not coincide when a is not an integer! Depending on the 
value of a, the result is similar to that seen in figure [7] on [71 This causes the 
function f(z) = not to be analytic on the whole plane; it is not analytic 
at the center of our disk, z = —r. To be technical again, we should omit 
the branch cut from the domain of the function and say that the function is 
only analytic on the region C \ (M~ U {0}) =C\{z : Re{z) < 0,lm(z) = 0}. 

However, another definition of an analytic function is a function which 
is continuously differentiable on a region G, and ^(-z"^) = az"~^, so clearly 
the function is differentiable. However, we must be careful with our choice of 
arguments (angles) to insure that the function is continuously differentiable. 
We will run into the same argument problem as in the definition of the 
logarithm, so we must make a branch cut, leading us to the same domain of 
analyticity mentioned in the previous paragraph. 

How can we adjust the definition of to create a function that possesses 
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Figure 7: The image of the disk from Figure [T] under the sequence of maps 
in Equation 13.11 on page [5] 

this continuous differentiabihty? We would like to arrange the branch cut 
so that when we follow the steps in Equation l3. li the branch cuts align after 
the exponentiation. Let us examine a formula for generating the values of 
and see how to address this problem. 

4 Computation and Formulae 

Computationally, what can we say about z"? How many distinct values in 
the complex plane should result? When q = we have a direct formula for 

the n*^ roots = z"' (see |2j, for example): 

zn = \z\ne^ " where 9 = arg(z), < k < n — 1 (4-1) 

which takes on n distinct values in the complex plane. So if we were to 
extend this to non-integer values for n, we would have an equation for a 
non-integral root of z: 

zv = \z\'^e^ ' where = arg(2;) (4-2) 

But how many values can k range over? Without loss of generality, let r/ = ^ 
with (p, g) = 1. Then we can interpret Equation 14.21 as 

zi = zi = \z\i ) where 9 = arg(z), Q < k < q — 1 (4.3) 

and we see that this formula will not cause the same value to appear until 
we have q distinct values for zi . This should make sense, since zi = {z'^)i . 
When we convert back to the rectangular coordinates, we see that these q 
points are exactly the same as the q points we compute using the g*^ roots 
of z^, but they appear in a different order. 
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This formula (Equation 14. 3p can be used to compute the a^^ powers 
(where a = |) of z directly: 

p p f (9 + 2fc7r)p \, 

Z9 = |z| 9ev 1 J where 9 = arg(2;), {) < k < q — 1 (4.4) 

However, this is not a one-to-one function; there are q distinct values. These 
values have arguments that cover more that an interval of 27r; in fact these 
q values have arguments that range over the interval [0, 2g7r)! This is an 
important note and we will return to it after an example. 
Example: Let us compute the | powers of 1 + i. We will number them 
ro, ri, . . . , r^. We have p = 5, (7 = 2, |1 + i| = ^/2, 6 = j. So we can use 
Equation 14.41 and we get: 

2 /(f+2ftvr)2V 

rfc=f\/2j%V J whereO<A:<4 

Now we evaluate this to get: 




Now compare this to the 5*^ roots of (1 + i)^ = 2i, which we will denote 
po, . . . , p4. We have \2i\ = 2, = ^, n = 5 so by Equation 14. 1^ we get: 



1 



+ 2fe7r 



Pk = \2\^e^ ^ ^ , where <k <A 
which evaluates to the following: 



k = ^ Pq = {2)h\ ' / = {2)h(^y 
k = l ^ pi = (2)te^ 




k = 3 ^ P3 = (2)56 



(5+8^) 



= 4 ^ P4 = (2)teV " J' = {2)h(^^ 
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We can see that there is a one-to-one correspondence between the values 
ro, . . . , r4 and the values po, . . . , p4^. 

po = (2)ie(T5)* = ro 

Pi= i2)h(^y =r3 

P2 = (2)5e'^ior = n 
P3 = [2)5 6^ 10 r = r4 

P4 = (2)5e(^)* =r2 

What we would like to develop is a way to avoid the branch cut problem 
seen in Figures [6] and [7] on page [71 We would like to arrange our domain so 
that the disks which cross the branch cut always realign (as in Figure 3] on 
pageE]). This will require a new domain for the function f(z) = z". 

5 A New Domain? 

For the function + c = e("'°s^) +c, with a = ^ a non-integral number, let 
us redefine the domain on a more general Riemann surface Wa- (Note that 
this set could also be called Wq, but we want to emphasize the dependence 
on a.) We construct the surface Wa by taking q copies of the plane, all 
slit along the negative real axis, each identified by the branch of arg(z) 
used on each sheet. The sheet we label as the 0*^ sheet has arguments in 
the interval (— 7r,7r], the sheet we label as the +1*^* sheet has arguments in 
(7r,37r], and so on up to the {q — sheet, which has arguments in the 
interval ((2g — 3)7r, (2(7 — l)7r]. In general the j*^ sheet has arguments in 
the interval ((2j — l)7r, {2j + l)7r] for any integer < j < q. The sheets are 
joined together at the negative real axis so that the argument is continuously 
rising as one travels in the positively oriented direction around the origin, 
and returning to the 0**^ sheet after the [q — 1)*** sheet. In this fashion, the 
surface is similar to the Riemann surface for z^/"^ with its finite number of 
sheets. An example is seen in Figure [8] on page [TOl 

We must examine the definitions of multiplication and addition in order 
to define F{z) = z"+c on this space properly. We must determine a method 
to keep track of the sheet we are on when we operate in this space. 

5.1 Multiplying and Adding in our Domain. 

Define (z,m) G (C,Z) as z = re^'^** = re^'^*(™+'^) where t G [m — i,m + |). 
This assumes that the branch cut is the negative real axis. We define mul- 
tiplication using the same idea as polar coordinates: given 

= rie^^^*^ = rie2^*("i+«i) and Z2 = rse^^^*^ = ^^e^^^^'^^+Ca) 



257r _ vr 

smce = — (mod ztt) 

10 2 ^ ^ 



337r _ ISvr 

smce = (mod Ztt) 

10 10 ^ ^ 
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Figure 8: The surface Wc 
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define 

{zi,mi){z2,m2) = {ziZ2,m^) = {rir2e^^^'^^^^^'^\m^) where ^1+^2 G 

Note that may or may not be equal to mi + 777-2; that depends on the 
values of ^1 and ^2- However, differs from mi +m2 by at most ±1. This 
allows us to define the non-integral powers of (z,m) (i.e., {z,m)°' = {z'^,m^,) 
for the appropriate value of m^,) in terms of the polar coordinates as well. In 
this context, keeping track of the sheet information is easy while multiplying. 

This means that the Riemann surface Wa completely removes the branch 
cut discontinuities at and 00 estabhshed in the definition of z". They are 
replaced by ramification points of Wa', i.e., points where is not locahy 
one-to-one. So we have established a multiplication that works when a is 
not an integer. 

However, when adding, keeping track of the sheet information is much 
harder. We want to define an addition operation (translation), called © to 
keep it straight, so that we can keep track of the sheet information our sum. 
So suppose that {z,mi) © (c, 7712) = {w,m^), and we define the values of w 
and as follows. Define w as the usual sum of 2; + c using rectangular 
coordinate addition. 

We use the following cases to determine m*; please note that QI is 
the first quadrant of the plane (i.e., the set {Re(2:) > 0,lm(z) > 0} ), 
QII is the second quadrant (the set {Re(z) < 0,lm(z) > 0}), QUI the 
third quadrant ({Re(2:) < 0,Im(2;) < 0}), and QIV is the fourth quadrant 
{{Ke{z) > 0, Im(z) < 0}). Since we want our operation to agree with the 
exponentiation operation, we want the addition to be continuous in 0, where 
z = re*^. This is referred to as "Counter-Clockwise Continuity" or CCC in 
[4]. While this algorithm seems hard to follow, the idea is simple: if we cross 
the negative half of the real axis then we change sheets; if we cross from QII 
to QUI, then we go up one sheet, but if we cross from QUI to QII, then we 
go down one sheet. 

Algorithm 5.1.1 To add the point (c, 0) to the point (z,m) in Wa, we 
choose the new sheet number m^, as follows: 

1. Re(c) = Xc = 0- There are three subcases: 

(a) Im(c) = Uc = 0. Then c = 0, w = z and m^ = m and we are 
done. 

(b) Im(c) = ?/c > 0. Then if z is in QI, QII or QIV, we have 

w = z + c in the sense of rectangular coordinate addition and 
m-t = m. However, if z is in QUI, we have two further subcases: 

i. I Im(c)| < I Im(2:)|. Then m^ = m 



771* 



1 1 
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ii. |Im(c)| > |Im(^;)|. Then m* = m — 1. Here is where we 
actually cross the branch cut. 

(c) Im(c) = yc < 0- Then if z is in QI, QUI or QIV, we have w = 
z + c in the rectangular coordinate sense and m* = m. However, 
if z is in QII, we have two further subcases: 

i. I Im(c)| < I Im(z)| . Then rritf = m. 

ii. |Im(c)| > |Im(2;)|. Then = m + 1. Here we cross the 
branch cut going the other way from the crossing above. 

2. Re(c) = Xc > 0. Again we have three subcases 

(a) Im(c) = ?/c = 0- Then m* = m for any z. 

(b) Im(c) = yc > 0. If z is in QI, QII, or QIV then m* = m, but if 
z is in QUI there are three subcases: 

i. If I Im(c)| < I Im(z)| then = m 

ii. If I Im(c)| > I Im(2;)| and \ Re(c)| > | Re(z)| then m* = m 
Hi. //|Im(c)| > |Im(2;)| and|Re(c)| < |Re(z)| then m^, = m — 1. 
This subcase, where c moves z up the imaginary axis, but not 
over far enough on the real axis to avoid the branch cut is the 
only one of these three subcases that crosses the branch cut. 

(c) Im(c) = i/c < 0. If z is in QI, QUI, or QIV, then m* = m, but if 

z is in QII there are three subcases: 

i. If I Im(c)| < I Im(z)|, then m,^ = m 

ii. If I Im(c)| > I Im(z)| and | Rc(c)| > | Re(z)| then rn^ = m 
Hi. //|Im(c)| > |Im(z)| and\Ke{c)\ < |Re(z)| thenm^=m + l. 

3. Re(c) = Xc < 0. Here also, there are three subcases: 

(a) Im(c) = yc = 0. Then m* = m for any z. 

(b) Iin(c) = t/c > 0. 

i. If z is in QI or QII, then m* = m. 

ii. If z is in QUI, there are two possibilities: Either |Im(c)| < 

I Im(z)| in which case m* =m, or\ Im(c)| > | Im(2;)| in which 
case nit, = m — 1 
Hi. If z is in QIV, there are three subcases 

A. If I Re(c)| < I Re(2;)| then = m 

B. If I Re(c)| > I Re(2;)| and \ Im(c)| < | Im(2;)| then m^ = m 

C. //|Re(c)| > |Re(z)| and |Iin(c)| > |Iin(2;)| then m* = 
m — 1 

( c) Im(c) = 2/c < 
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i. If z is in QUI or QIV then m.^ = m 

a. If z is in QII the there are two possibilities: Either \ Im(c)| < 
|Im(2;)| in which case = m, or |Im(c)| > |Im(z)| then 
m=K = m + 1 
Hi. if z is in QI, then there are three subcases: 

A. if I Re(c)| < I Re(2:)| then = m 

B. if I Re(c)| > I Re(2;)| and \ Im(c)| < | Im(2;)| then = m 

C. i/|Re(c)| > |Re(z)| and |Im(c)| > |Im(2;)| then = 
m + 1 

If we want to generalize to the case where we are adding and 
{z2,m2), the process is nearly identical. We wish to preserve the idea from 
addition of vectors that the sum of two vectors lies "between" the vectors 
(i.e., the parallelogram law). The major change in the method above is that 
m* is now based on [I^itHMJ ^nd then adjusted by ±1 depending on the 
cases above. The only main concern is that the point is on every sheet. 
Hence \i Z1+Z2 = then the sheet is irrelevant. By default, the sheet should 
be left as [^Hiii^J ^ in order to simplify the definition. This now allows us 
to add any two complex numbers in the Riemann surface setting. We will 
define neighborhoods of zero topologically: any ball of radius e around zero, 
covering all sheets, is a neighborhood of zero. 

The major change is that addition is now a noncommutative operation! 
It is order dependent, as we see in the following example. 

Example 1: Suppose we are working on a 4-sheeted space, like W15/4. 
Let (zi,mi) = (— 2 + i,2) and {z2,m2) = (—1 — 3i,2). Converting these 
to polar coordinates, we get zi ^ ^e°-^^24i67ri ^ ^g0.6024i67rj^ 

So when we compute [zi + Z2,?n.*) we get 2:1+22 = —3 — 2i and m^, = 3 
(by case 3(c)(ii) above). However, when we compute (22 + zi,m^,) we get 
22 + 2i = — 3 — 2« and = 2 (by case 3(b) (ii) above). □ 

We run into a big problem: this operation is not continuous. 

Counter-example 1: Let us consider an e -neighborhood of the point 

2o = (^,1), under the map 2 1— )• 2 + i, on a three sheeted surface like Ws. 
^ 3 

Then when we translate, using Algorithm 15. l.H the points in Arn{Re(2) < 0} 
are shifted to sheet 0, while the points in n {Re(2) > 0} remain on sheet 
1. Hence is "sheared" into two half-disks: one containing the point (^, 1), 
and one containing the point (^,0) in its closure. □ 

This algorithm is specifically constructed to match the choice of the neg- 
ative real axis as the branch cut. This is not the only way to construct the 
surface Wa and still maintain continuity in the argument, t, in our construc- 
tion. We could use a branch cut along the positive real axis, or along any 
curve between the points and 00, just as any of these curves define a valid 
branch of the logarithm. But this will not change the continuity problem. 
In fact, any rigid translation on Wa is inherently discontinuous. Once we 
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decide that our domain is a ramified Riemann surface like Wa, translation 
becomes discontinuous at the ramification point (0 in our case). 

Theorem 5.1.2 Let C be the branch cut from to oo of the surface Wa- 
For any c G Wa, there exists a point zq and a neighborhood N = N{zq) 
such that the set N + c = {z + c:z£ N} is not a connected set (using a 
translation defined similarly to Alaorithm \5.1.1\ with C as the branch cut). 
In particular, we may choose zq such that G {zq + tc, < t < 1}. 

Proof: Fix a and c. Choose zq such that G {zq + tc, < t < 1}. Then 
the ray from the origin through c divides = N{zq) into two subsets Nl 
and Njj, each of which is a half-disk. When is translated, points in one 
of the half-disks change sheets, while points in the other half-disk do not 
change sheets. Regardless of which sheet the point zq is mapped onto, half 
of will be separated from zq. If we approach zq in A'^^^, then for the sake 
of continuity, zq should be mapped onto the sheet containing A'^^;,; but the 
same argument holds true for Nfj. But the point zq is itself only mapped 
onto one sheet. Hence the translation map is not continuous. □ 

This answers a question raised by Robert Corless in his E.C.C.A.D. pre- 
sentation [3]: "Can a Riemann surface variable be coded? What will the 
operations be on it?" (this was also discussed in detail in the Appendix A to 
E. Kaltofen's paper [lOj-) Unfortunately, we are answering "No, a Riemann 
surface variable cannot be coded. There is no continuous addition operation 
on such a variable." 

We might also choose a different algorithm for choosing m*. For example, 
we can define m* based on the sign of Im(c): 

Algorithm 5.1.3 To add the point (c, to the point (z,mi) in Wa, we 
choose the new sheet number as follows: 

1. //Im(c) > then = mi + m2 + 1. 

2. //Im(c) = then m* = mi + m2. 

3. //Im(c) < then = mi + m2 — 1. 

But this algorithm for choosing sheets is not continuous either. 

Counter-example 2: Fix a = | and c = (i,0), so Wa has two sheets, 
with branch cut along the negative real axis. Consider an e-neighborhood 
A^ of the point (—1,0). Now half of N is on sheet and half of A^ is on 
sheet 1. Hence when we add i using Algorithm I5.1.3| the neighborhood A^ 
is "sheared" into two parts, since the points on sheet 1 are moved to sheet 
and vice-versa. A^ will become two half-disks centered about the points 
(i, 0) and {i, 1). 

By defining the function have two points in C where 

F is not conformal, namely and oo. This is caused by the branching of 
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log(z) and the failure of the branch cuts to hue back up after multiplying by 
a (recall the Figures |6] and [7]) . Any translation by c cannot prevent these 
points from staying ramified. Hence ramification at zero implies either c is 
also a ramification point or ^ i— )• z + c is not continuous. Since ramification at 
c requires infinitely many ramification points (since c is chosen arbitrarily), 
this isn't a good plan. 

6 Looking Forward 

We come to the following negative result: the operation of addition is in- 
compatible with the operation of exponentiation for non-integer exponents. 
So in order to study the fractals generated by functions of this form, we 
must redefine them on the plane and make sure that the definitions match 
the definitions for the integer case as closely as possible. We have examined 
these sets in a separate article (|13|). 

Many of the characterizations of fractals for polynomials will fail to carry 
over to these functions and that the discontinuity we observed in Section [2] 
becomes more important than ever. This discontinuity is observed when an 
algorithm for generating fractals is given non-integer exponents. Figure [9l on 
pagell6| is an example of the Julia set for the function f{z) = z^'^ -|- \i. This 
image was generated using the Fractint program (version 19.6) found at the 
web site given in [5j (using the fractal type julzpower), and the calculations 
were carried out using only the principal branch of the logarithm. The 
parameters for the julzpower fractal type are: the real and imaginary parts 
of the parameter (0 -|- 0.5i); the real and imaginary parts of the exponent 
(2.5 -|- Oi); the bailout test and the bailout value, both of which are set to 
for the defaults (modulus for the test and 4 for the value). 
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Figure 9: The Julia set for the function f{z) 
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